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ABSTRACT 

We present results from a study to determine whether relations - established in the local Universe - between 
the mass of supermassive black holes (SMBHs) and their host galaxies are in place at higher redshifts. We 
identify a well-constructed sample of 18 X-ray- selected, broad-line Active Galactic Nuclei (AGN) in the Ex- 
tended Chandra Deep Field South - Survey with 0.5 < z < 1.2. This redshift range is chosen to ensure that 
HST imaging is available with at least two filters that bracket the 4000 A break thus providing reliable stel- 
lar mass estimates of the host galaxy by accounting for both young and old stellar populations. We compute 
single-epoch, virial black hole masses from optical spectra using the broad Mgll emission line. For essentially 
all galaxies in our sample, their total stellar mass content agrees remarkably well, given their BH masses, with 
local relations of inactive galaxies and active SMBHs. We further decompose the total stellar mass into bulge 
and disk components separately with full knowledge of the HST point- spread-function. We find that ~ 80% of 
the sample is consistent with the local M B h _ ^*, Bulge relation even with 72% of the host galaxies showing the 
presence of a disk. In particular, bulge dominated hosts are more aligned with the local relation than those with 
prominent disks. We further discuss the possible physical mechanisms that are capable building up the stellar 
mass of the bulge from an extended disk of stars over the subsequent eight Gyrs. 
Subject headings: galaxies: evolution — galaxies: active 



1. INTRODUCTION 

A determination of the physical mechanisms through which 
supermassive black holes are built up at the centers of galax- 
ies have been one of the ke y issues in astrophysics (see 
Kormend y & Richstondll995l) . Such processes are thought 
to further provide a link black hole growth and the formation 
of the bulges of their host galaxies based on both observa- 
tions and theory. Correlations betwee n the mass of the central 
black hole and absolute magnitude (Magorr ian et al.l 119981: 
iMarconi & Huntl2003l: fe aring & Rix 2004]), and/or stellar ve- 
locity dispersion (Gebhard t et al.l 120001; iMerritt & Ferraresel 
1200 1|) of the spheroidal component indicate that the mass 
ratio between a SMBH and its bulge is constant over a 
wide dynamic range in mass (e.g. Mg///Mg M /^= 0.0014; 
lHaring & Rixl (12004 : hereafter M BH - M Bu i ge relation). We 
will refer to this relation as the local relation. 

Over the past years, several studies have addressed whether 
there is an evolution of the mass relations between the central 
black hole and its host galaxy. Such studies must rely upon 
galaxies with accreting SMBHs (i.e., Active Galactic Nuclei; 
AGN) since the region of influence surrounding black holes 
cannot be resolved at higher redshifts. While those hidden by 
obscuration (i.e., type 2 AGNs) give a rather clean view of 
their host galaxy, unobscured (type 1) AGN are the only sys- 
tems for which black hole masses can be measured. Although, 
an estimate of the mass of the host bulge is challenging due 
to the glare of a luminous AGN that only gets more difficult 
at high redshift. Fortunately, optical imaging from space with 
HST can be used to d isentangle the light between an AGN and 
its host galaxy (e.g ISanchez et al. 2004; Jah nke et al.ll2004l 
120091: iBennert et aLll2011bt iCisternas et all 120 111) due to the 
high spatial resolution and well understood point spread func- 
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tion. Alternatively, it is also possible to measure the stellar ve- 
locity dispersion from optical spectra for less luminous AGNs 
(IWoo etal.ll2008h : this method requires high signal-to-noise 
spectra that limits its application to high redshift AGNs. 

Even if the host galaxy is resolved only limited spec- 
tral coverage is usually available to estimate stellar masses. 
Single-band studies are therefore restricted to the black hole 
mass - luminosity relation or have to make assu mptions on 
the m as s-to-light ratio of th e host galaxy (see iPeng et al.l 
I2006allbl : iDecarli et al.ll2010allbh . iMerloni et all (|2010|) imple- 
mented a new approach to measure the stellar mass content of 
AGN host galaxies through template fitting of the broad-band 
photometric sp ectral energy distribution (Brusaetal. 2009; 
IXue etal.ll20T0h . With this approach. IMerloni et al.l (120101) 
estimate the total stellar mass content which provides only an 
upper limit to the bulge mass. Bennert et al. (201 1) take a sig- 
nificant step forward by using the multi-ba nd HST data avail- 
able in the GOODS (iGiavalisco et al.ll2004l) fields to decom- 
pose the AGN and host galaxy light including a bulge com- 
ponent tractable through multiple filter bandpasses. Unfortu- 
nately, the sample is selected to be a redshifts (1 < z < 2) for 
which the optical imaging falls below the rest-frame 4000 A 
break. Surprisingly, the aforementioned studies find elevated 
black hole masses as compared to eith er the bul ge component 
(Woo et al. 2008; B ennert et al. 120 11 bl) or total (Merloni et al 



1201 Ol) stella r mass of their host gal axy. Recently, Jahnk e et al 



(|2009|) and ICisternas et al.l (|2011|) report that the mass ratio 
between the black hole and the total stellar mass of its host 
galaxy is similar to local values possibly indication of an un- 
dermassive bulge. 

Even with the considerable effort achieved to date, there 
are several challenges that need to be met in order to accu- 
rately determine the evolution of the Mbh - ^*, Bulge at higher 
redshift. First, the decomposition of optical light is more dif- 
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ficult due to the strong surface brightness dimming of the host 
galaxy as compared to the AGN. To mitigate this effect, high 
resolution imaging with high signal-to- noise is needed to ad- 
equately resolve the host galaxy especially for bright AGN. 
Equally important, at least, one rest- frame optical color and a 
luminosity is n eeded to constra in the stellar mass content of 
the host galaxy dBell et al.l2003l) . A color that covers the 4000 
A break provides a good estimator on the underlying stellar 
mass-to-light ratio. It is worth highlighting that the 4000 A 
break moves out of the optical filter bands at z > 1.2 thus re- 
quiring deep high-resolution NIR imaging. Furthermore, due 
to the limited physical resolution at high redshift and the fact 
that galaxies become more compact, it may be challenging to 
classify galaxies morphologically such as distinguishing be- 
tween disturbed and undisturbed hosts. 

A determination of the M B h~^*, Bulge relation using AGN 
samples, also requires an assessment of the possible bi- 
ases originating from sel ection of AGN (see Salviander et al.1 
120071: lLauer et ail |2007|) . While quiescent galaxies are se- 
lected by their magnitude or luminosity, active galaxies (e.g., 
unobscured, broad-line AGN) are often selected by their op- 
tical nuclear luminosity or magnitude. The bias introduced 
by the luminosity (i.e., mass) limit will have a stronger ef- 
fect at the high mass end of the black hole mass func- 
tion which is strongly decreasing. Offsets from the local 
M B h - *, Bulge relat ion seen in samples of luminous AGN with 
massive BHs (e.g . jMerloni et al.ll2010l;lBennert et al Jl201 lbh 
Peng et al. 2006a : b) may be explained by such a bias. There- 
fore, a sample selected at lower luminosities that fall well be- 
low the knee of the black hole mass function should be less 
impacted by such a bias. 

In this study, we determine the Mbh -^*, Total and Mbh~ 
Af*, Bulge relations at 0.5 < z < 1.2 using a sample of 18 
X-ray selected broad-line AGN (BLAGN) from the Ex- 
tended Chandra Deep Field - South Su rvey. B ased 
on HST/ ACS imaging from G EMS dRix et al J [2004b and 
GOODS ( Giavalisco et al. 2004), we measure the stellar mass 
content of their host galaxies including the bulge component. 
We specifically focus on this redshift range so that there is 
at least one HST band above and below the 4000 A break 
thus providing a rest- frame color required for accurate con- 
version of light to mass. Black hole masses are determined 
using single epoch virial mass estimation based on the Mgll 
emission line. In Section 2, we describe our sample. In Sec- 
tion 3, we describe our analysis of the HST/ACS data that 
involves the image decomposition of the total light into AGN, 
bulge and disk decomposition, and stellar mass estimation. 
Black hole masses are fully detailed in Section 4. Section 
5 and 6 presents the results including a discussion of the re- 
lations between the mass of the SMBH and their total/bulge 
stellar mass. Finally, in Section 7 we give a summary of the 
results. Throughout this paper we assume a flat cosmology 
with #o = 70 kms" 1 Mpc" 1 , Q M = 0.3 and ^ A = 0.7. 

2. AGN SAMPLE 

Currently, broad-line (type 1) AGNs provide the only 
means to establish the relation between BH mass and galaxy 
mass beyond the local universe. This is due to the fact 
that BH mass measurements rely upon a determination of 
the velocity widths of gas in the vicin ity of the BH as 
provided by broad emission lines (e.g., Kaspi et al.1 12005 
Vestergaard & Peterson 2006). High-resolution imaging (best 
if taken from space) can then be used to detect the extended 
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Fig. 1. — Optical R-bmd magnitude versus broad-band X-ray flux (0.5-8 
keV; units of ergs cm 2 s _1 ) for all X-ray sources from Lehmer et al. 2005 
(grey circles). All sources classified as BLAGN have been marked with a 
black dot. Our final sample selection is shown as magenta solid squares. 
Objects falling below the typic al relation for type I AGNs are the obscured 
(type II) AGNs (see Figure 5 of Silverman et al. 20®. 



emission from the underlying host galaxy. There have been 
numerous studies of the host galaxies of type 1 AGNs using 
such techniques (e.g JJahnke et a l. 2004, 2009; Sanche z"et al.l 
2004: lBennert et alJ^OlfblilCisternas et al.ll2011l) . 

We aim to take advantage of type 1 AGNs that are found 
in X-ray surveys, such as the Chandra Deep Field South 
- Survey, that reach faint depths. These X-ray sources are 
likely to have a wide range in their optical properties that 
includes those of lower luminosity both missed in optically- 
selected samples such as SDSS, and more favorable for the 
study of their host galaxy. There are numerous papers on the 
host galaxies of X-ra y selected AGN that may be of inter- 
est to the reader (e.g.. iGrogin et al.ll2005t [Pierce et al.ll2QQ7t 
lAmmons et al . 2009; Silve rman et al.ll2008l) . 

Here, w e specifically se l ect typ e 1 AGNs from the com- 
pilation of ISilverman et al.l (2010) that provide spectroscopic 
redshifts and classification of X-ray sources i n the the Ex- 
tende d Chandra Deep Field South - Survey (|Lehmer et al.l 
120051) . These are objects with at least one broad emission 
line having a FWHM greater than 2000 km s -1 . We further 
require that an available spectrum has a good enough qual- 
ity to perform our emission line fitting procedure to estimate 
virial black hole masses. 

We then demand that each type 1 AGN has been observed 
by H ST. The ECDF S is covered by the G EMS dRix et al.1 
120041) and GOODS dGiavalisco et al.l 120041) surveys in the 
central area. GEMS consists of imaging in two optical 
HST filters (ACS F606W and F850LP) while GOODS has 
four filters (ACS F435W, F606W, F775W and F850LP). Un- 
fortunately, some sources are located on the outskirts of 
the ECDFS and therefore no HST coverage is available. 
Even though, extensive ground based data is available of the 
full ECDFS are a from various observi ng campaigns (e.g., 
MUSYC survey; Cardamone et al. 20lo|). we choose to avoid 
any biases that may appear due to the inclusion low resolu- 
tion data. We further stress that is essential to have at least 
two filters that bracket the 4000 A break in the rest- frame of 
the host galaxy for accurate estimation of the mass-to-light ra- 
tio and the stellar mass (see the following section). To do so, 
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we elect to restrict the type 1 AGN sample to0.5<z<1.2 
that allows us to determine accurate rest- frame B-V colors for 
the entire sample. In addition, we apply the same selection to 
the deeper 2Msec catalog (Luo et al. 2008) and identify one 
additional source. Our final sample consists of 18 type 1 AGN 
with half falling in the GOODS area and the other half within 
the GEMS field. In Figure 1, we show the distribution of X- 
ray flux and R-band optical magnitude of Chandra sources 
and highlight those within our type 1 AGN sample. It is ap- 
parent that the sample spans about two dex in both X-ray flux 
or luminosity, and optical brightness. The final sample cov- 
ers the full region of the fx—R plane as the overall BLAGN 
sample. 

3. OBSERVED HOST GALAXY PROPERTIES 
3.1. AGN -host decomposition and bulge correction 

The first step to obtain information on the host galaxy 
is to remove the contribution of the AGN component from 
the broad-band HST images. The separation of galaxy light 
from that of the nuclear point source in luminous AGN is 
challenging even at lower redshifts since the AGN can out- 
shine the host galaxy by several magnitudes. Objects in our 
study have lower luminosities (due to their X-ray selection 
with deep observations) thus the contamination from the point 
source is substantially weaker compared to similar studies us- 
ing optically- seleted quasars. 

Knowledge of the point spread function (PSF) at the posi- 
tion of the AGN is crucial for our further analysis. The ACS 
PSF in the GEMS survey is known to vary across the field 
( Jahnke et al. 2004). Therefore, we create a local PSF for 
each AGN by averaging all stars within a radius of 60 arcsec 
around the target. Each high S/N PSF consists of about 30-40 
stars. The remaining uncertainties between individual stars 
are included in the variance frame as an additional contribu- 
tion from the rms image of each PSF. In the GOODS fields, 
the estimation of a proper PSF is more difficult. Each tile has 
only a limited number of unsaturated stars (5- 20) with strong 
spatial variations, in some cases, between stars located in the 
center and at the edges of each tile. For most of our objects, 
we used the same strategy as for GEMS by creating a local 
mean PSF for each AGN position. In three cases, the AGN 
was close to the edge of the tile thus we used the nearest star 
as our PSF reference. 

We use GALFIT dPeng etalJl200l [20loh to fit the two- 
dimensional light distribution of each AGN with a point 
source model represented by an empirical PSF plus a Sersic 
model for the host galaxy. The nucleus component is either an 
average PSF created from various field stars around the target, 
or a single PSF star as described above. The decomposition 
of the images is done in several steps. First, we conserva- 
tively subtract a PSF scaled to the flux contained in a small 
aperture (typically 2 pixels) around the central pixel. For the 
second step, we perform a full decomposition by adding a Ser- 
sic model as a second component. We then optimize the fit to 
minimize the residuals. This step requires several iterations 
of the fit using different starting values to ensure convergence 
to a global minimum in the parameter space. If necessary, we 
add further components to fit asymmetries in the host (e.g., 
arm structures). Neighboring galaxies are either masked out 
or fitted simultaneously (see ID-333 for such an example in 
the bottom panel of Figure 2) to avoid flux spilling over from 
one object to the other. If the host galaxy flux is below 2%, 
we decide that the host galaxy is unresolved. 



In Figure [2 we show three representative examples of 
the image decomposition procedure for objects with differ- 
ent nuclear-to-host (N/H £| ratio and bulge-to-total (B/T) stel- 
lar light ratio. As shown, these cases demonstrate the effec- 
tiveness of both the short (F606W) and long (F850LP) wave- 
length HST imaging. The top panel shows ID- 158 that ex- 
hibits a nuclear component attributed to the AGN and a clearly 
extended component characterized by a Sersic index of 4.2. 
With no discernable disk, the morphology is determined to be 
that of an early-type galaxy. In the middle panels, the host 
galaxy of AGN (ID-417) is well resolved above our detection 
limit even though it has a high N/H ratio. The Sersic index is 
found to be n=2.1 that does not favor either a simple early or 
late type morphology. Only through a decomposition of the 
bulge and disk components (as described below) can we de- 
termine whether this object is truly bulge or disk dominated. 
The third example (ID-333) shown in the bottom panels has 
a galaxy contribution with a Sersic index of 1.25 indicating a 
strong disk contribution to the overall morphology. In Table 
[U we list the sample properties and the results of our fitting 
routine. In Figure [3] we show the PSF- subtracted host images 
of the entire sample. 

We estimate the uncertainties of our measurements through 
a series of simulations. We create artificial AGN images using 
empirical PSFs and host galaxy models superimposed with 
artificial noise to match the flux levels measured in the real 
images. We estimate statistical errors on the host galaxy and 
nuclear magnitude by comparing the input and output values 
for our fit parameters. Host galaxy apparent magnitude, ra- 
dius and morphology are extracted from the Sersic model fits. 
Since the Sersic index can be underestimated, in some cases, 
using GALFIT, especially when the nuclear-to-host (N /H) ra- 
tio is high (see lSanchez et al J 120041: [Kim et alJl2008albh . we 
use the simulations to correct the Sersic index. Here, we also 
want to point out the importance of the HST data again spe- 
cially in cases such as ID-333 (See Fig. 2) where the AGN is 
strongly contaminated by a nearby companion that is hardly 
resolved in ground-based data. 

A direct bulge/disk decomposition is only possible for ob- 
jects which have low N/H ratios (typically N/H < 2). We 
find 5/18 host galaxies to have n sers i CiCorrecte d > 3; therefore, 
we classify them as truly bulge dominated. If the single Sersic 
fit of the host galaxy indicates the possible presence of a disk 
component with n sers i C ^ correcte d < 3 (13/18 objects), we refit the 
galaxy with two Sersic models each representing the disk and 
the bulge. We put limits on the Sersic Index (0.5 < n < 1.5 
for the disk and 3 < n < 5 for the bulge) of each component to 
achieve an effectively reduced chi square of the residuals as 
compared to using single values typical for a disk (n=l) and 
bulge (n=4) components. In the end, we are able to directly 
decompose 14/18 host galaxies in the F850LP filter into either 
a purely bulge or bulge+disk component. Some bulge+disk 
component fits failed due to the disturbed morphology of the 
host galaxy (i.e. ID-271 or ID-516) even though the AGN 
was weak (N/H < 1). Since the F850LP filter provides the 
best contrast between nuclear and host component, we can use 
the best fit parameters (i.e. disk/bulge radius, position angle, 
axis ratios) as constraint in the bluer filter bands with typically 
higher N/H ratios. The best fit radii for the bulge components 
range from 0.6-6 kpc and are consistent with the typical sizes 

2 The nuclear-to-host ratio (N/H) is simply the flux attributed to the AGN 
(N) divided by that of the host galaxy (H). Both determined through decom- 
position of the HST images. 
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Fig. 2. — Three examples of HST image decomposition and surface profile fitting. Top panels: AGN ID- 158 having a bulge-to-total (B/T) light ratio equal to 
1.0 based on the F606W filter band. The four images are as follows: (i) the original image (upper left), (ii) the host galaxy after removing the point source (upper 
right), (iii) the best fit model (lower left), and (iv) the residual after subtraction of the best fit model from the original image. The scaling is the same in all images. 
On the right (upper panel) we show the surface brightness profiles of the various components (i.e. the original data - filled gray circle, the nucleus model/PSF 
red dotted line, the host galaxy after removal of the PSF (green pentagon), the best fit single sersic model - dashed line and the overall fit as a solid line) In the 
lower panel we show the residual after subtraction of the best fit model profile from the data. The single Sersic fit indicates a early type galaxy; therefore, no 
further decomposition into bulge and disk components is necessary. Middle panel: AGN ID-417 (B/T=0.3) in the ACS/F850LP filter band. Bottom panel: 
AGN ID-333 (B/T=0.35) in ACS/F850LP filter band. Due to the large contribution of the companion to the surface brightness profile, we model the companion 
separately and show its component and contribution to the overall surface brightness profile. Both ID-417 and ID-333 have host galaxy emission that can be 
decomposed into bulge and disk components. 
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TABLE 1 

Sample: Target List and results from the analysis. The different columns show (1) Object ID taken from the Lehmer et al., (2), (3) 

RA DEC COORDINATES, (4) REDSHIFT, (5) ABSOLUTE V-BBAND MAGNITUDE, (6) REST-FRAME U-V COLOR OF THE HOST GALAXY ,(7) SERSIC INDEX 
FROM SINGLE SERSIC FIT,(8) HALF-LIGHT RADIUS OF THE BULGE AND DISK IN ARCSEC INF850LP (9) TOTAL STELLAR MASS, (10) BULGE-TO-TOTAL 
LUMINOSITY RATIO, (11) FWHM OF THE MGlI EMISSION LINE, (12) CONTINUUM LUMINOSITY AT 3000A, (13) BH MASS, (14) EDDINGTON RATIO, 
(15) NUCLEAR-TO-HOST RATIO IN F850LP, (16) SURVEY FIELD: GEMS [1] & GOODS [2] 



Source ID RA Dec z M bost U-V Sersic r bu ige/rdisk M tot [M ] B/T FWHM C L 30 oo log Mbh e N/H Field 

[1] [2] [3] [4] [5] [6] Index [7] [8] [9] [10] [11] " [12] [13] [14] [15] [16] 



158 


52.942418 


-27.952407 


0.717 


-21.70 b 


0.78 


4.0 


0.83/- 


10.66 


1.00 


3.24 


43.22 


7.08 


0.07 


0.05 


[1] 


170 


52.949839 


-27.845910 


1.065 


-20.33 T 


0.04 


2.2 


0.07/0.12 


9.80 


0.50 


3.12 


42.85 


6.98 


0.05 


0.11 


[1] 


250 


53.001530 


-27.722073 


1.037 


-21.20^ 


0.84 


7.7 


0.13/- 


10.19 


1.00 


4.12 


44.01 


7.64 


0.09 


1.00 


[2] 


271^ 


53.125255 


-27.756535 


0.960 


-21.08" 


0.76 


1.4 


-/- 


10.49 


0.23 


2.72 


43.86 


7.26 


0.19 


0.54 


[2] 


273" 


53.016951 


-27.623708 


0.970 


-20.39 b 


0.96 


4.3 


0.22/- 


10.29 


1.00 


6.96 


43.84 


8.13 


0.03 


1.22 


[1] 


305 


53.036121 


-27.792822 


0.544 


-22.47 T 


0.58 


4.0 


0.77/- 


11.01 


1.00 


5.80 


44.98 


8.52 


0.14 


1.51 


[2] 


333 


53.057789 


-27.602162 


1.044 


-21.06 T 


0.41 


1.1 


0.08/0.14 


10.33 


0.35 


6.40 


43.29 


7.80 


0.02 


0.19 


[1] 


339 


53.062421 


-27.857514 


0.675 


-22.00 T 


1.04 


5.5 


0.32/- 


10.96 


1.00 


9.81 


42.60 


7.74 


0.00 


<0.03 


[2] 


348 


53.071454 


-27.717531 


0.569 


-21.5 F 


0.72 


1.2 


0.21/0.24 


10.55 


0.40 


5.82 


43.91 


7.78 


0.04 


0.08 


[2] 


375 


53.110383 


-27.676530 


1.031 


-22.4 l" 


-0.12 


1.1 


-/- 


10.55 


0.10 


3.02 


45.18 


7.92 


0.63 


3.14 


[2] 


379 


53.112521 


-27.684732 


0.737 


-22.49 


0.02 


2.4 


-/- 


10.66 


0.32 


10.30 


45.05 


9.04 


0.05 


3.28 


[2] 


413 


53.156053 


-27.666706 


0.664 


-20.40 b 


0.81 


1.9 


0.30/0.36 


10.19 


0.41 


2.29 


43.38 


6.95 


0.16 


0.29 


[2] 


417 


53.158807 


-27.662460 


0.837 


-20.36 T 


1.02 


2.1 


0.10/0.37 


10.39 


0.33 


5.22 


44.67 


8.25 


0.12 


5.50 


[2] 


465 


53.199639 


-27.696625 


0.740 


-21.603 


0.87 


2.2 


0.40/0.51 


10.76 


0.24 


5.60 


43.78 


7.92 


0.04 


0.22 


[1] 


516 


53.246074 


-27.727665 


0.733 


-22.48" 


0.46 


2.3 


-/- 


10.88 


0.67 


6.08 


43.44 


7.80 


0.02 


0.08 


[1] 


540 


53.256378 


-27.761801 


0.622 


-22.03 b 


0.79 


1.6 


0.32/0.67 


10.88 


0.49 


5.18 


43.07 


7.52 


0.02 


<0.03 


[1] 


597 


53.292483 


-27.811635 


1.034 


-21.93 T 


0.92 


2.8 


0.27/0.66 


10.91 


0.77 


7.22 


43.54 


8.01 


0.02 


0.18 


[1] 


712 


53.370529 


-27.944765 


0.841 


-22.98^ 


0.96 


2.3 


0.21/0.96 


11.34 


0.70 


6.44 
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a [l] ID taken from Giacconi et al. 2002 

b [2] direct bulge estimate through either B or B+D fit to imaging data 
c [3] units of 1000 km/s 
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FIG. 3. — Host galaxy images after removing the point source using GAL- 
FIT in the F850LP filter band for the whole sample. We have marked objects 
for which we could get a direct bulge estimate through either a pure bulge (B) 
or bulge+disk fit (B+D) 



of elliptical galaxies at similar redshifts ( Truj il loet al .1 12007 ) . 
In four cases, the radius of the bulge component is less than 
3.5 pixel and the fit can be treated as an upper limit. But we 
also want to caution that the radii are more sensitive as a free 
parameter of the fit than the fluxes of the components. 

We can further check whether the inclusion of a pseudo- 
bulge component provides a reasonable fit to the surface 
brightness profile. To do so, we allow the bulge component 
to have a lower Sersic index. While it is challenging to dis- 
tinguish between a pseudo-bulge and a classical bulge with 
the data in hand, we do find stronger residuals in the nuclear 
region if we fit with a pseudo-bulge component. These new 
fits lead to only small changes in the fluxes of each compo- 
nent (<0.1 mag) and bulge-to-disk ratio (<13%), hence only 
a small impact on the stellar mass estimates. The poorer 
fits, seen when incorporating a pseudo-bulge component, may 
indicate that SMBHs are more directly related to classical 
bulges than pseudo-bulges (Kormendy et al. 20TT|). 

For N/H>3, we estimate a bulge correction through exten- 
sive simulations rather than a direct measurement. For each 
filter band starting with the longest wavelength due to lower 
contrast between AGN and host galaxy, we create a set of ar- 
tificial AGN+host images using the best fit parameters for the 
nucleus component. The host galaxy consists of two Sersic 
models, one for the disk and one for the bulge component. 
We vary the free parameters to mimic a broad range of bulges 
and disks. We add appropriate noise measured from the HST 
images and fit the artificial images with a single Sersic model 
plus PSF model. We compare the fit solutions of the simula- 
tions with our single Sersic fit of the real data and select the 
B/T ratios of all model fits that recover the original fit parame- 
ters within their uncertainties. We fit all filter bands separately 
to account for possible color gradients between the disk and 
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Fig. 4.— B/T distribution of ID-417 (top) and ID-375 (bottom) extracted 
from our simulations matching the properties of the host galaxy using a single 
Sersic component. The insets shows the host galaxy in F850LP after removal 
of the nuclear point source. 



the bulge. Due to higher N/H ratios at shorter wavelength and 
typically lower S/N (fainter objects) we constrain some of the 
free parameters such as size, ellipticity and centroid position 
to the range of solutions found in the F850LP filter band. In 
Figure |U we show the B/T distributions recovered for two of 
our host galaxies (ID-417 and ID-375) together with an image 
of the host galaxy. Figure [4] also shows the good agreement 
between the recovered B/T ratio and the host morphology. For 
example ID-375 has a low B/T ratio of about 0.1 indicating 
the presence of a dominating disk and the image of the host 
galaxy shows some spiral structures with very little light con- 
centration in the center of the galaxy. The mean value of the 
B/T distribution and its uncertainty is then used to estimate 
the bulge mass. The error bars on the photometry are typi- 
cally larger by 0.15-0.3 mag. 

3.2. Total/Bulge Stellar Mass Estimates 

Our main goal is the estimation of the stellar mass content 
of each host galaxy and its bulge component in the sample 
by converting the rest frame optical colors into mass-to-light 
ratios. For targets in the GEMS area, we have only a sin- 
gle optical color while for GOODS we have multiple colors 
based on four filter bands. This method has been success- 
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fully employed in several studies on AGN host galaxies (e.g., 
ISchramm et al.ll2QQ8h IJahnke et al.ll2QQ4ll2QQ9l:lsinchez et al.1 
12004 . As demonstrated by Bell et al. (2003) for a variety of 
star formation histories, the stellar mass-to-light ratio (M/L) 
can be robustly predicted from the B-V color. We adopt their 
formula based on a Chabrier IMF: 



log 10 (Af/L y ) = -0.728 + (1.305 x (B-V)) , 



(1) 



where M/Ly is given in solar units. The choice of an IMF 
has a systematic effect on the final mass estimation. Using a 
Salpeter IMF would typically increase our mass-to-light ratios 
by factor of 1.4. 

The ECDFS area is well covered by broad-band photome- 
try from various instruments ranging from the ultraviolet to 
the infrared. The available photometry provides another ap- 
proach to estimate the total stellar mass content of the host 
galaxies through direct SED fitting (as shown by Merloni et al. 
2010). Although this would be a powerful alternative, the data 
suffers from additional uncertainties such as source confusion 
in the ground based data or variability of the sources due to 
multi-epoch data. In any case, we decided to implement SED 
fitting only as a consistency check on our total stellar mass es- 
timates based on the HST data. For the procedure, we use our 
own algorithm based on a Levenberg-Marquar t y 2 minimiza- 
tion. We use a set of SED model templates (Maraston 2005) 
with declining star formation histories based on a Kroupa 
IMF (which gives similar results as a Chabri er IMF), solar 
metall icity and a dust extinction law following ICalzetti et al.l 
(2000). We make use of the broad-band image decomposi- 
tion (AGN+host) based on the HST results to constrain the 
template models that includes the photometric errors in each 
filter band. First, we fit a template AGN model (Richards et 
al. 2006) to the photometry of the nucleus obtained from the 
decomposition and subtract this from the total (ground- and/or 
space-based) photometry. Next, we fit the residual fluxes with 
either a single template or a two component template. Strong 
contamination from unresolved sources in the ground based 
data (i.e. in ID-333) are taken into account and subtracted sep- 
arately using the HST photometry as an additional constraint 
for the companion template model. To estimate errors on the 
stellar mass, we use a Monte Carlo approach. We vary the 
observed flux in each bandpass by a random number which 
is Gaussian distributed with a sigma defined by the flux er- 
ror. We generate 100 simulated SEDs and recompute the fit. 
Masses from both approaches typically agree within 0.13 dex. 
In Figure 6, we show four examples (ID-339, ID- 170, ID-465 
and ID-250) of the SED decomposition compared to the sin- 
gle epoch spectrum used for the BH mass estimation. These 
four AGN represent different level of AGN and host stellar 
continuum throughout our sample. In addition, these objects 
also have different B/T ratios: 1.0 (ID-339, ID-250), 0.5 (ID- 
170) and 0.24 (ID-465). These examples illustrate the clear 
advantage gained from having the HST photometry. Only 
with HST resolution, we can constrain the flux of the AGN 
and the host galaxy including the bulg e and disk c ompo nents 
separately. The mass estimates from IB ell et al.1 ([2003) and 
SED fitting typically agree within 0.15 dex. 

In Figure [51 we plot the U-V rest- frame color versus the to- 
tal stellar mass. Nearly all hosts concentrate at the high mass 
end and below the red sequence (i.e., the green valley). The 
masses of our host galaxies are comparable to typical red se- 
quence galaxies but the colors of the host indicate a popula- 
tion of recently formed stars. 



10 1 



10 1 



data 

■ nucleus model 
host model 

■ composite model 




4000 



5000 6000 7000 8000 

observed wavelength [ 



9000 



10 1 



10 1 



■ data 

nucleus model 

host model 

~_ composite model 


i i i i 

ID-465 
z=0.740 






"Til \s 

kiN i i i 



4000 



5000 6000 7000 8000 

observed wavelength [ 



9000 



10" 1 



10 1 



10 1 



data 

■ nucleus model 
host model 

■ composite model 




4000 5000 6000 7000 8000 9000 

observed wavelength [ 



data 

■ nucleus model 
host model 

■ composite model 



ID-339 
z=0.675 




4000 



5000 6000 7000 8000 

observed wavelength [ 



9000 



FIG. 5. — Examples of the SED template fitting for different nuclear-to- 
host ratios (ID-339JD-250; ID-170; ID-465). Black squares show the to- 
tal photometry measured in the bands covered by HST (either F606W and 
F850LP for GEMS or F435W, F606W, F775W, F850LP for GOODS). Red 
and blue squares show the results from the image decomposition in the HST 
filter bands representing the nuclear component (red) and the host component 
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models with the magenta solid line showing their sum. 
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4. BLACK HOLE MASSES 

We measure black hole masses for our entire type 1 AGN 
sample using single-epoch spectra that provide both a veloc- 
ity width of a broad emission line and the monochromatic lu- 
minosity of the continuum. We use optical sp ectra acquired 
mainly from the followup of X-ray sources (|Szokoly et al.l 
120041: [Silverman et"aT1l2Q 10). We supplement these with spec- 
tra taken with FORS2 on the VLT but not yet publicly avail- 
able. 

Several prescriptions to estimate black hole masses 
are available from the literature using vario us emis- 
sion lines such as H/3, M^II o r CIV (jKaspi et all 



2000; Vestergaar d & Petersonl llM ICollin et a l. ~ 2006 



McLure & Jarv is 2002). Due to the redshift range of our sam- 



ple and optical spectroscopic coverage, we use the Mgll emis- 
sion line to estimate virial black hole masses in all cases. Al- 
though most of the black hole mass calibrations are based 
on reverberation mapping data of Hf3 several studies have 
shown that there is good agreement between the mass esti- 
mates based on Mgll and the B aimer lines (H(3, Ha) out to 
high redshifts (Shen & Liu 2012; Matsuoka et al. 2012) by 
combining optical and NIR spectroscopy. The prescription 
for es timating black hole mass as given in iMcLure & Jarvisl 
(2002) is implemented; although, we recognize that similar 
recipes are available elsewhere (Kong et al. 2006, McGill et 
al. 2008, Wang et al. 2009) with each of these agreeing es- 
sentially to within 0.2-0.3 dex. 

We perform an iterative least- squares minimization to fit the 
Mgll line for each AGN to measure its line width. Our proce- 
dure is a modified version of the one used in Gavignaud et al.l 
(2008). The number of components to fit the line depends 
on the characteristics of the objects and quality of the data. 
We fit the region around the emission line using a model that 
includes a pseudo-continuum and one or two Gaussian com- 
ponents to characterize the line profile. We find that for the lo- 
cal continuum a powerl aw+broadened Fe-template (provided 
by M. Vestergaard: see Vestergaard & Wilkes 2001) gives the 
best results. Specially the strength of the Fe-emission in the 
wings of the Mgll line can vary strongly (see ID-250 for 
strong Fe and ID- 170 for very weak Fe) and affects the out- 
come of the fit. We try to both minimize the number of model 



components and optimize the residuals around the emission 
line. We either interpolate over absorption features or mask 
them out. A FWHM of the line profile is determined using ei- 
ther a one or two component Gaussian model. We have teste d 
the same algorithm on the sample from lMerloni et al.l (|2010|) . 
Even tough we find some scatter for the individual fits, there 
is no systematic offset in the final black hole mass estimates. 
Two of our objects overlap with the study from lBennert et al.l 
(|2011b|) : our mass estimates agree within 0.1 dex using the 
same recipe. 

In the next step, we measure the continuum luminosity at 
3000 A required to estimate a radius to the BLR. For lumi- 
nous AGN (Lboi > 45) the continuum luminosity can be di- 
rectly measured from the spectrum due to the typically low 
impact of the host galaxy. For our sample, we find that in sev- 
eral cases, there is a significant host galaxy contribution that 
must be taken into account (see Fig. [5]). Therefore, we de- 
cided to measure the monochromatic luminosity at 3000 A by 
decomposing the HST/ACS images. The procedure enables 
us to isolate the AGN (i.e., nuclear) emission from its host 
galaxy m ost effectively. We th en fit an average quasar SED 
template (Richards et al.ll2006|) . accounting for dust attenua- 
tion to estimate the intrinsic continuum luminosity at 3000 A. 
We find that the continuum luminosity based on HST imag- 
ing agrees with that determined from the decomposition of 
the broad-band SED to within 5%. Monte Carlo realizations 
using the uncertainties of the FWHM and L3000 measurements 
enable us to estimate the uncertainties on the black hole mass 
in addition to 0.4 dex uncertainty inherent in the scaling rela- 
tions. In Figured we present examples of the fits to the broad 
emission lines in six AGN with different quality of data. A 
summary of the results of our line fits are shown in Table 1 . 

5. RESULTS 

5.1. The BH Mass- Total Stellar Mass Relation 

We first present the relation between black hole mass (Mbh) 
and total stellar mass (M* Total) in Figure [8] (/e/Y panel). From 
the distribution of data points, it is apparent that our sample 
does not have the dynamic range in either stellar mass or black 
hole mass to establish both a slope and normalization simul- 
taneously of a linear fit. Fortunately, we can compare with 
the local relation establishe d using inactive galax ies (mainly 
ellipticals or SO) as done by Haring & Rix| (12004)) and deter- 
mine whether an offset exists. We find that essentially all of 
our AGN fall along the local Mbh Bulge relation. It is im- 
portant to highlight that the bulge mass is equivalent to the to- 
tal stellar mass for the local comparison sample. To be more 
specific, we find that 17/18 objects, considering their la er- 
rors, are consistent with the typical region of 0.3 dex scat- 
ter around the best fi t local relation having a slope of 1.12 
( Haring & Rix 2004). Given our limitations in mass cover- 
age as mentioned above, we fit a linear regression model to 
our data while fixing the slope to the value given above thus 
determining only the normalization. We fi nd the best-fit nor - 
malization to be 8.31 by using FITEXY (iPress et al.lll993b . 
which estimates the parameters of a linear fit while consider- 
ing errors on both variables. The fit is affected by the single 
target offset from the relation. If excluded for no obvious rea- 
son, the constant would be 8.24. With a simple Monte Carlo 
test, we can reject the null hypothesis that the two samples are 
significantly different. While the local inactive sample is es- 
tablished using dynamical masses, we do not expect these to 
differ substantially from the stellar masses; this is in fact the 
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Fig. 7. — Multicomponent fit to broad Mgll emission line for six representative objects from our sample. The top panel shows the spectral range around the 
emission line. The best fit model is indicated as a red solid line. The different components are as indicated: Fe-emission (blue), pseudo-continuum (purple), 
Gaussian components (green and brown). The residual (data- fit) is shown in the lower panel. In the upper right corner of each panel, we show the black hole 
mass distribution computed from our Monte Carlo tests based on the uncertainties of the line width and continuum luminosity measurements. The objects are 
(from top left to bottom right): ID-170,ID-250,ID-413,ID-417,ID-339, and ID-379 



case as demonstrated in Ben nert et al.l (1201 lab . 

Ideally, we would like to compare our sample with a local 
sample of active SMBH s with stellar mass me asurements of 
their hosts. The work of Benner t et al.l (|201 la) allows such a 
direct comparison. We show these data in Figure [8]as marked 
by small black circles. Carrying out the same fit as for the 
ECDFS AGNs, we find the best fit constant to be 8.30 for 
the local AGNs. We use the t otal stellar ma ss for the regres- 
sion fit of both active samples ( Benne rt et al.l (1201 1 ah . ECDFS 
AGNs) and find no significant deviation between them in the 
M B u -M^Totai relation. 

Our result agrees well with findings of recent studies of 
the Mbh -M* .Total relation at high redshift. In particular, 
JahnkeetaLl J2009) use a similar technique of decomposing 
HST images and converting rest- frame optical colors into stel- 



lar mass-to-light ratios based on a sample of AGNs at z > 1 
in COSMOS with NICMOS coverage. Their sample consists 
of ten objects with seven for which they achieve a decom- 
position in multiple bands and find no offset in black hole 
mass given their total stellar masses. Our study effectively im- 
proves the statistics by a factor of 2.5 and fills in a gap in red- 
shift coverage (see Figure [8] right pa nel). In addition, these 
results are supported by the findings of ICisternas et al.l (1201 1|) 
who explored the same relation on a sample of BLAGN at 
0.3 < z < 0.9 from the COSMOS survey; although, only one 
HST band is available to constrain the stellar mass content of 
the host galaxy. 

Taken toge t her, these studies ( Jahn keeit al.l [2009; 
ICisternas et al.l 1201 1|) . including our own, clearly con- 
trast with other works at high redshift that claim an increasing 
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offset in black hole mass for a given stellar mass. In the 
right panel of Figure [5] we show the redshift evolution of 
M B h-^*, Total ratio compared to various other studies probing 
the same relation. Some studies are using different mass 
estimators fo r the black hole m asses or stellar masses. 
For example , Merloni et al.l (|2010|) use the prescription of 
iMcGill et all (120081) for their black hole mass estimations 
and assume a Salpeter IMF for their stellar mass estimates. 
When necessary, we convert the masses of different studies to 
the pr escription based on the formula from iMcLure & Jarvisl 
(2004) and a Chabrier IMF for the stellar mass estimates. 
In case of iMerloni et al.l (|2010), the corrections have only 
a marginal effect on the M B h - ^*,Totai relation. Based 
on our results, we cannot confirm or rule out a stronger 
evolution at higher redshift (z > 1.5). In particular, our mean 
bolometric lumin osity is log Lpqi =4 4.7 while the higher 
redshift sample of Mer loni et al.l (1201 0|) is at log Lt i = 45.5. 
As a consequence, the mean BH mass is shifted to higher 
masses and therefore a direct comparison with these objects 
and any trend implied by the data might be biased by the 
differences in the sample properties. It is worth highlighting 
that our results are likely to be less biased due to selection 
since our BH masses are typically below 1O 9 M , the knee in 
the black hole mass function; we may be effective ly avoiding 
the problems fully presented in lLauer et al.l (l2007|) . 

5.2. The BH Mass-Bulge Stellar Mass Relation 

While the total stel lar mass is well-de t ermined us- 
ing different methods ([Schramm et all 120081: iJahnke et al.l 
2009; Me rloni et alJ2010tlCisternas et al.l20TTtlBennert et al.l 
201 ltj), we usually do not know how much of the total mass is 
present in the bulge. As stated above, only 5/18 of our AGN 
hosts have a Sersic Index n > 3 indicating a purely bulge- 
dominated host galaxy. We make the assumption that for 
these objects the total mass is the same as the bulge mass. For 
the remainder, we estimate the bulge contribution to the total 
mass by corrections to the total mass by accounting for the 
contribution of the disk. Applying the same cut at n < 3, we 
find that ~72% of the host galaxies show a disk component. 
Although the fraction is in good agre ement with the results 
presented by ISchawinski et al.l (|2011|) on a sample of X-ray 
selected AGN in the Chandra Deep Field South at 2 < z < 3. 
Although, we draw a different conclusion on the importance 
of the disk component, in terms of the mass contribution to the 
total mass. Our bulge/disk decomposition shows that, even 
though a disk is present, the mass of the central bulge can 
still dominate the total mass of the host galaxy. The different 
redshift regimes might play an important role since there is 
ab out 3-5 Gyr of galaxy ev olution between our study and that 
of ISchawinski et al.l (|2011|) . Using the B/T ratio to divide our 
sample into bulge and disk dominated systems, we find that 
~ 50% of the sample has a significant bulge component with 
B/T > 0.5; this can even be true for objects with a surface 
brightness profile of the host galaxy described by a fit with a 
Sersic index of ~2. 

We can now establish the M B h - ^*, Bulge relation at 0.5 < 
z < 1.2. In Figure |9l we plot the Mbh -M* : Bulge relation 
and compare our results with the sample of in active galaxies 
from H aring & Rix (2004) and local AGN from Bennert et al. 
(|2011a|) . The stellar mass measurements for the local AGN 
allow a more direct comparison with our sample than the dy- 
namical masses of Haring & Rix (2004). We find that the 
mass distributions for all three samples are very similar with 
each other. This can be clearly seen in a histogram of the mass 



ratio (log M B H/M* iBu ige) shown in the top panel of Figure [TOl 
where there is no significant difference in the median value. 
Overall, we find that 78% of the AGNs are consistent with the 
local relation. If we consider the single object undergoing a 
clear major merger (ID-333), there are only three objects that 
are significantly offset from the local relation. If we artifi- 
cially move this object onto the relation, then 83% of AGNs 
in our sample are consistent with the local relation. We inter- 
pret this as evidence for a black hole-bulge relation, at these 
redshifts, to be similar to the local relation. Interestingly, we 
do find additional scatter in our sample compared to that in the 
local distributions. We further note that there are no objects 
well below the M B h - ^*,Buige relation. 

We can further investigate where high-z AGN lie in respect 
to the local relation as a function of their bulge-to-total ratio. 
All objects with B/T > 0.5 fall nicely onto the local relation 
(see Figure [9] and the bottom panel of [Tot . They are also 
the most massive objects in the sample in terms of their bulge 
mass. Objects with a B/T < 0.5 are clearly separated in bulge 
mass (from bulge-dominated objects) and the majority are still 
in good agreement with the local relation. Only four objects 
have under massive bulges considering their 1 a error bars in- 
cluding ID-333 which has a massive companion that might 
move the whole system onto the relation after the merger. 

6. DISCUSSION 

An important question for SMBHs and their host galax- 
ies is their subsequent evolution in the black hole - bulge 
mass plane. As previously mentioned, 83% of the bulges 
in our sample are already massive enough that their M B h - 
Bulge ratio agrees well with that seen in inactive galax- 
ies today (see Figure [9]). We illustrate this further in Fig- 
ure [TOl (top panel), by comparing the distribution of the 
M B h -M*^ Bulge ratio between various samples. Interestingly, 
there are some outliers with undermassive bulges, relative to 
their BH mass, that are preferentially disk dominated galax- 
ies. In the bottom panel of Figure [TUl we compare the distri- 
butions of this ratio for the bulge and disk dominated subsam- 
ples separately to the distribution of the local AGNs. Even 
though the number statistics are small, we find no difference 
for the bulge dominated subsample by looking at their me- 
dian ratios. The situation is different for objects in the disk- 
dominated subsample. While some objects overlap with the 
distribution M B n —M^ Bulge ratios of the local AGN, the me- 
dian ratio of the disk dominated subsample is shifted by 0.5 
dex towards a higher ratio. When comparing bulges of simi- 
lar mass (log Bulge < 10.5), local AGN host galaxies have 
a smaller offset (~0.25 dex). On the other hand for the same 
mass matched subsample which includes 22/25 objects in the 
local AGN sample and 12/18 from our sample, we find that 
the local AGN sample contains only ~30% disk dominated 
systems while our subsample contains ~ 80% disk dominated 
systems. Within the AGN population, we may be witnessing 
both a migration onto the local relation and a morphological 
transformation with cosmic time. We recognize that selec- 
tion effects may impact such comparisons. Ideally, we want 
to have an AGN sample spanning a wide baseline in redshift 
with equivalent selection, BH mass indicators and sufficient 
statistics. 

This leads to the question how these host galaxies can grow 
their stellar bulge mass to match the bulge masses seen to- 
day. One possible track could be the event of a major merger 
that leads ultimately to a significant increase in stellar bulge 
mass. Mergers are seen to play a role in black hole growth for 
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Fig. 8. — Left: Mbh _ ^*,Totai relation for our sample of intermediate redshift AGN. The sample has been divided into bulge (red squares) and disk dominated 
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and the local AGN (stellar mass) sample from Bennert et al. (201 la). The dashed line represents the best fit from Haring & Rix 2004 with a 0.3 dex scatter shown 
as the grey shaded area. Right: Redshift evolution of the relation in comparision with several other studies taken from the literature. The dashed line shows the 
mean constant ratio from the local relation by Haring & Rix (2004). 
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Fig. 9. — Mbh -A^*,Buige relation for our sample of intermediate redshift AGN. Symbols are the same as in Figure[8] 



similar X-ray selected samples (Silv erman et al.ll20TT)) . Out 
of our 18 AGN, only one (ID-333) shows signs of an ongoing 
major merger. Even though other host galaxies do show some 
signs of minor merger activity, we conclude that the growth 
of the bulge through a major merger event in the near future 
is not certain. On the other hand, the good agreement of the 
AGN host galaxy Mbh-M* 5 Total relation with the local relation 
clearly shows that all the mass needed to put our host galax- 
ies onto the local Mbh ~M* ; Bulge relation is already in plac e 
within these galaxies at redshift z ~ 1 (Jahn ke et al.l 12009). 
Therefore, mass transfer from the disk to the bulge is necces- 
sary to grow their bulges. Any bulge growth through internal 



processes has to overcome the mass growth of the black hole 
otherwise the galaxy would just move on a diagonal track in 
the M B h - M*,Bui g e relation. 

While the BHs in their active phase are growing, we can 
also investigate how the host galaxy is growing in stellar 
mass by looking at their individual growth rates (i.e., SFR) 
and compare these to the BH growth rates. We estimate 
star-formation rates based on the UV continuum from our 
best fit SED models and converted these into growth rates 
(SFR/M s teii). In FigureQj] we compare the growth rates of the 
host galaxies with the growth rates of the BHs as determined 
by M/M B H- M is determined from L^/ = cMc 2 . To estimate 
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log M BH /M st 




log M BH /M stell bul 

Fig. 10. — Top panel: Mbh -Af*, Bulge ratio distribution for our sample of 
intermediate redshift AGN (red solid) compared to the sampl e of inactive 
(black dotted) and active (green dashed) galaxies from Bennert et al. (201 la) 
using s tellar host galaxy masses. The inactive sample from Bennert et al. 
( 2011b) is based on the group 1 sample from Marconi & Hunt (2003) which 
overlaps also with the sample of Haring & Rix 2004. Bottom panel: Com- 
parision of the Mbh —M* ,Buige ratio distribution for the bulge dominated 
(Lb/Lt > 0.5) subsample shown as red solid line, the disk dominated 
(Lb/Lt < 0.5) subsample shown as blu e solid line. The dat a is superim- 
posed on the sample of local AGN from Bennert et al. (2011a) showing the 
sample separation into bulge (red dashed) and disk (blue dashed) dominated 
systems. 



the bolometric luminosities L^ / and Edding ton ratios, we 
use th e luminosity dependent corrections from Hopkins et al. 
( 2007) applied to our derived continuum luminosities at 3000 
A. We find that apparently the BHs gain mass much stronger 
than the host galaxies by a factor of ~30. These relative 
growt h rates are broadly consistent with th at seen in obscured 
AGN dNetzeri2009HSirverman et al.l2009l) . Su ch an offset im- 
plies that the typical duty cycle of an AGN (see Martini ( 2004) 
for an overview) during which it can grow its BH mass effi- 
ciently must be short enough (typically 10 7 - 10 8 yr) to pre- 
vent a significant vertical movement in the black hole mass 
- bulge mass plane. If the growth rates are extrapolated over 
a period of lGyr, the host galaxies do not gain much stel- 
lar mass from the present level of star formation. As pre- 
viously mentioned, only one object (ID-333) shows a possi- 
ble major merger due to the presence of a more massive but 
inactive companion. While some objects show signs of mi- 
nor merging activity (i.e. ID-712JD-271), and we cannot ex- 
clude that we miss further minor merger events due to their 




log Growth Rate [yr 

Fig. 11. — Comparision of the growth rate distributions of the host galaxy 
(blue; M^^ Tota i /SFR) estimated from the UV-SFR and the BHs (red; Mbh x 
(dM/dty 1 ) of our sample 



low surface brightness, the stellar mass gain is expected to 
be low. Assuming the current growth rates and ignoring a 
possible major merger, all galaxies except one would need 
more than >~ lGyr to move more than 0.3 dex in the Mbh- 
Total- Therefore, we do not expect much evolution over the 
next lGyr for the majority of our sample. 

7. SUMMARY 

We have performed a detailed analysis of a sample of 18 
type 1 AGN host galaxies at 0.5 < z < 1.2 to estimate their 
stellar mass content and explore the relation between the mass 
of the central BH and the mass of the host galaxy. Our sam- 
ple is of moderate-luminosity due to a selection based ini- 
tially on their X-ray emission as detected with the Extended 
Chandra Deep Field - South Survey. This results in a sample 
having black hole masses below the knee of the black hole 
mass function thus mitigating biases (Lau er et aT1l2007|) seen 
in other samples to date. For the chosen redshift range, HST 
imaging is available with at least two filters that bracket the 
4000 A break thus providing reliable stellar mass estimates of 
the host galaxy by accounting for both young and old stellar 
populations. We have estimated bulge masses for all galax- 
ies through either direct decomposition of the imaging data 
into a bulge or bulge plus disk component, or through simu- 
lations where artificial host galaxies with different B/T ratios 
are compared to single Sersic fits of the host galaxy. We are 
now able to look separately into their relation of the BH mass 
with either total stellar mass content or bulge mass after the 
contribution from the disk is removed. 

We find that the relation between M B h and Total is in very 
good agreement with the local Mbh -M^ Bulge which has been 
reported by several studies so far. From our morphological 
analysis and decomposition of bulge and disk components, 
we can quantify the fraction of bulge dominated objects with 
B/T > 0.5 to be 50% while 72% of the sample shows the 
presence of a disk component which is a significantly higher 
fraction than for a stellar mass matched local AGN sample. 
Even though the bulge mass is shifted towards lower masses 
given their BH mass in some cases, we find that ~80% of the 
sample is in agreement with the local Mbh and M* Bulge rela- 
tion given their la error bars. We further compare the growth 
rates of the host galaxy and their BHs and find that assuming 
the present SFR and accretion rates (while ignoring possible 
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major merger events), only one AGN in our sample would 
move more than 0.3 dex over the next lGyr. We highlight 
that bulge dominated galaxies are well in place at z ~ 1 on the 
local Mbh - ^*, Bulge relation. There is a significant fraction 
(20%) of our sample that is disk dominated and above the lo- 
cal relation which is not seen in either local inactive or active 
galaxy samples. For these galaxies to grow their bulges and 
align themselves on the local relation a physical mechanism 
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